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Introduction 

LanHEP has been designed as part of the CompHEP package [1], a software for automatic 
calculations in high energy physics. CompHEP allows symbolic computation of the matrix 
element squared of any process with up to 6 incoming and outgoing particles for a given physical 
model {i.e. a model defined by a set of Feynman rules as a table of vertices in the momentum 
representation) and then numerical calculation of cross-sections and various distributions. 

The main purpose of the new option given by LanHEP is to easily introduce new models. 
Although this job is rather straightforward and can be done manually, it requires careful cal- 
culations and in modern theories with many particles and vertices, such as supersymmetric 
models, can lead to errors and misprints. 

LanHEP makes possible the generation of Feynman rules for propagators and vertices start- 
ing from the Lagrangian defined by a user in some simple format very similar to canonical 
coordinate representation. It is possible to use multiplets and substitution rules, for exam- 
ple for covariant derivative. Supersymmetric theories can be described using superpotential 
formalism and 2-component fermion notation. 

The user should prepare a text file with description of all Lagrangian terms in the format 
close to the form used in standard publications. The output of the program is Feynman rules 
in CompHEP model tables, or in FeynArts [2] format, or LaTeX table. 

The LanHEP software is written in C programming language. The first version [4] was 
released in August 1996. 

1 Getting started with LanHEP 
1.1 QED 

We start with a simple exercise, illustrating the main ideas and features of LanHEP. The first 
physical model is Quantum Electrodynamics. 



model QED/1. 




parameter ee=0. 31333: 'elementary electric 


charge ' . 


spinor el/El : (electron, mass me=0. 000511) 




vector A/A: (photon) . 




let F"mu'"nu=deriv"mu*A''nu-deriv~nu*A~mu . 




Iterm -l/4*(F''mu''nu)**2 - l/2*(deriv''mu*A' 


'mu) **2 . 


Iterm El*(i*gamma*deriv+me)*el . 




Iterm ee*El*gamma*A*el . 





Figure 1: LanHEP input file for the generation of QED Feynman rules 



The QED Lagrangian is 

jOqed = -^F^.F'^" + eY{id^ + geA^)e - mee 
and the gauge fixing term in Feynman gauge has the form 

Here e{x) is the spinor electron-positron field, m is the electron mass, A^^x) is the vector photon 
field, Pi^^ — d^A" — d^A^, and is the elementary electric charge. 



2 



The LanHEP input file to generate the Feynman rules for QED is shown in Fig. 1. 

First of all, the input file consists of statements. Each statement begins with one of the 
reserved keywords and ends by a full-stop '.' symbol. 

First line says that this is a model with the name QED and number 1. This information is sup- 
plied for CompHEP, the name QED will be displayed in its list of models. In CompHEP package 
each model is described by four files: 'varsN.mdl', 'funcN.mdl', 'prtclsN.mdl', 'IgrngnN.mdl' , 
where N is the very number specified in the model statement. 

The model statement stands first in the input file. If this statement is absent, LanHEP 
does not generate the four standard CompHEP files, but just builds the model and prints a 
diagnostic, if errors are found. 

The second line in the input file contains declaration of the model parameter, denoting 
elementary electric charge Qe as ee. For each parameter used in the model one should declare 
its numeric value and an optional comment (it is also used in CompHEP menus). 

The next two lines declare particles. Statement names spinor, vector correspond to the 
particle spin. So, we declare electron denoted by el (the corresponding antiparticle name is 
El) and photon denoted by A (with antiparticle name being A, since the antiparticle for photon 
is identical to particle). 

After the particle name we give in brackets some options. The first one is the full name of 
the particle, used in CompHEP; the second option declares the mass of this particle. 

The let statement in the next line declares the substitution rule for symbol F. 

Predefined name deriv, which is reserved for the derivative will be replaced after the 
Fourier transformation by the momentum of the particle multiplied by —i. 

The rest of the lines describe terms in the Lagrangian. Here the reserved name gamma denotes 
Dirac's 7-matrices. 

One can see that the indices are written separated with the caret symbol Note that in 
the last two lines we have omitted indices. It means that LanHEP restores omitted indices 
automatically. Really, one can type the last term in the full format: 

Iterm ee*El"a*gamma~a"b"mu*A"mu*el''b. 

It corresponds to geCaJab^b^n with all indices written. Note that the order of objects in the 
monomial is important to restore indices automatically. 



1.2 QCD 

Now let us consider the case of Quantum Chromodynamics. The Lagrangian for the gluon 
fields reads 

Lym = -\f'''''F;,, 



where 



^ ii.v — (^fi^v ~ '^u'^fj. ~ 9sJ '^^'^u^ 



G^(a;) is the gluon field, gg is the strong coupling constant and Z"**^ are purely imaginary 
structure constants of the SU{3) color group. 

The quark kinetic term and its interaction with the gluon has the form 



Lf = qa^d^qi + gsXtjqa'^qjG''^: 



where Af, are Gell-Mann matrices. 
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Gauge fixing terms in Feynman gauge togetlier witli the corresponding Faddev- Popov ghost 
term are 

where (c, c) are unphysical ghost fields. 

The corresponding LanHEP input file is shown in Fig. 2. 



model QCD/2. 

parameter gg=l . 117: 'Strong coupling'. 

spinor q/Q: (quark, mass mq=0.01, color c3) . 

vector G/G:(gluon, color c8, gauge). 

let F~mu"nu~a = deriv"nu*G~mu~a - deriv"mu*G"nu"a - 

gg*f _SU3''a''b''c*G''mu~b*G~nu''c . 
Iterm -F**2/4-(deriv*G) **2/2 . 
Iterm Q*(i*gamma*deriv+mq)*q. 
Iterm i*gg*f_SU3*ccghost(G)*G*deriv*ghost(G) . 
Iterm gg*Q*gamma*lainbda*G*q. 

Figure 2: Input file for the generation of QCD Feynman rules 



Table 1: QCD Feynman rules generated by LanHEP in LaTeX output format 



Fields in the vertex 


Variational derivative of Lagrangian by fields 


Gfxp Vq V r 

Qap Qbq G 

Gpp Gjyq Gpj, 

n c c c 

^ HP ^fq ^ pr ^as 


~9sPsfpqr 

gsfpqMg''" p'2g^'' - f^g"' + vig^"" + v'ig'"' - p\g^') 

gl{g^'g'"'fpqifrst - g'^'g'^'fpqtfrst + g^'g'^'fprtfqst 
+g^''g"'fpstfqrt - g^^g^'fprtfqst - g^'g'^'fpstfqrt) 



Since QCD uses objects with color indices, one has to declare the indices of these objects. 
There are three types of color indices supported by LanHEP. These types are referred as color 
c3 (color triplets), color c3b (color antitriplets) , and color c8 (color octets). One can see 
that color c3 index type appears among the options in the quark q declaration, and the 
color c8 one in the gluon G declaration. Antiquark Q has got color index of type color c3b 
as antiparticle to quark. LanHEP allows contraction of an index of type color c3 only with 
another index of type color c3b, and two indices of type color c8. Of course, in Lagrangian 
terms each index has to be contracted with its partner, since the Lagrangian has to be a scalar. 

LanHEP allows also to use in the Lagrangian terms a predefined symbol lambda with 
the three indices of types color c3, color c3b, color c8 corresponding to Gcll-Mann A- 
matrices. Symbol f_SU3 denotes the structure constant /"^"^ of color SU{3) group (all three 
indices have the type color c8). 

Option gauge in the declaration of G allows to use names ghost (G) and ccghost(G) for the 
ghost fields c and c in Lagrangian terms and in let statements. 

Table 1 shows Feynman rules generated by LanHEP in LaTeX format after processing the 
input file presented in Fig. 2. Four gluon vertex rule is indicated in the last fine. Note that the 
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output in CompHEP format has no 4-gluon vertex explicitly; it is expressed effectively through 
3- leg vertices by a constant propagator of some auxiliary field (see section 10 for more details). 



1.3 Higgs sector of the Standard Model 



model Higgs/1 . 

parameter EE = 0.31333 : 'Electromagnetic coupling constant', 
SW = 0.4740 : 'sin of the Weinberg angle (PDG-94) ' , 
CW = Sqrt(l-SW**2) : 'cos of the Weinberg angle'. 

let g=EE/SW, gl=EE/CW. 

vector A/A: (photon, gauge), 

Z/Z:('Z boson', mass MZ = 91.187, gauge), 
'W+'/'W-': CW boson', mass MW = MZ*CW, gauge). 

scalar H/H: (Higgs, mass MH = 200, width wH = 1.461). 

let B = -SW*Z+CW*A. 

let W = {'W+', CW*Z+SW*A, 'W-'}. 

let phi = { -i*gsb('W+') , (vev(2*MW/EE*SW)+H+i*gsb(Z) )/Sqrt2 }, 
Phi = anti (phi) . 

Iterm -2*lambda*(phi*aiiti(phi)-v**2/2)**2 where 

lambda= (g*MH/MW) **2/ 16 , v=2*MW*SW/EE . 
let D''a''b~mu =(deriv~mu+i*gl/2*B"mu)*delta(2)''a''b 

+i*g/2*taupm"a"b"c*W"mu~c , 
Dc~a~b~mu=(deriv~mu-i*gl/2*B~mu)*delta(2) "a"b 

-i*g/2*taupm"a"b"c*anti (W) "mu"c . 
Iterm D''a''b''mu*phi"b*Dc''a"c"mu*Phi"c . 



The third example illustrates using the multiplcts in the framework of LanHEP. Let us 
consider the Higgs sector of the Standard Model. The Higgs doublet can be defined as 



where is the sinus of the weak angle, H is the Higgs field, Zf and Wj^ are goldstone bosons 
corresponding to Z and gauge fields. The Higgs potential reads as 



where A = {gMff/MwY/W, and the vaccuum expectation value v — 2Mw/g (here and below 

g = e/syj, g' = e/c^.) 

The gauge interaction of a Higgs doublet is given by the term (D^$)(£'^$)*, where 



Figure 3: Input file for the Higgs sector of the Standard Model 




D^^d^ + ig'Bj2 + igfW, 



Here B is C/(l) singlet and W is SU{2) triplet, 




= c^Z^, + s^A^ 



\ 



I 
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Table 2: LanHEP output: Feynman rules for Higgs self- interaction 



Fields in the vertex 


Variational derivative or Lagrangian by nelds 


ff H H 

11 11 11 


3 eMH"- 
2 MwSu} 


H Wt 


1 eMH'^ 

2 M\ySw 


n //IP / 1 I? 


1 eMH^ 

2 M^^Sw 


H H H H 


3 e^MH'^ 

4 Mw-2s„2 


12 n vv p vv p 


1 e'HlH^ 

4 A%-2su,2 


H H Zp Zp 


1 e^MH^ 
4 Mw'''sw'^ 


Wp Wp 


1 e^MH^ 

2 Mw^s^^ 


wt Wp Zp Zp 


1 e^MH^ 
4 Mw^s™^ 


Zp Zp Zp Zp 


3 e^MH'^ 

4 Mw^s,,,'^ 



and f is the vector (r"^, t^,t~). 

The above model can be represented by the LanHEP code shown in Table 3. 

New features in this example include the definition of special symbols for coupling constants 
g, g', gauge and Higgs fields, and the covariaiit derivative by means of let statement. Note 
that in the declaration for the fields we have used dummy indices (vector and isospin). The 
multiplets are defined by the components in the curly brackets. 

Note that the option gauge in the declaration of gauge fields allows to use the name gsb(Z) 
and gsb('W+') for the goldstone bosons. 

2 Structure of LanHEP input file 

The LanHEP input file is the sequence of statements, each starts with a special identifier (such 
as parameter, Iterm etc) and ends with the full-stop '.' symbol. Statement can occupy several 
lines in the input file. 

This section is aimed to clarify the syntax of LanHEP input files, i.e. the structure of the 
statements. 

2.1 Constants and identifiers 

First of all, each word in any statement is either an identifier or a constant. 

Indentfiers are the names of particles, parameters etc. Examples of identifiers from the 
previous section arc particle names 

el El A q Q G 

The first word in each statement is also an identifier, defining the function which this statement 
performs. The identifiers are usually combinations of letters and digits starting with a letter. 
If an identifier does not respect this rule, it should be quoted. For example, the names of W"^ 
bosons must be written as 'W+' and 'W-', since they contain '+' and '-' symbols. 
Constants can be classified as 
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Table 3: LanHEP output: Feynman rules for Higgs gauge interaction 



Fields in the vertex 


Variational derivative of Lagrangian by fields 


A W+ 


ltiIViy[r 


A W+ W~ 
vyp yy V 




■^n yy F yy F 


^\P2 P3) 


n yy yy v 


eMw 


n yy yy f 




£1 vv F yy 


2 s^. yP2 Pi ) 


H Z Z 




H Z,, Zp 


I " (Pi Ps) 


l^+„ Z„ 


r y 


W+. Wp Zp 


2 s-w ^■^•^ ^ ' 


w+ w-^ z. 


ieMwSw 


Wt W-^ Zp 


\tiP^l-P^^) 


W+ Wp Z^ 


_i(l_2^( 


A, Wp 




A^ H Wp 




A^ H W+ 


2s^y 


Wp Zp 


2s.ui^ 


A^ W+ Zp 


_ i ei f.nv 


A^ Wp Z, 


(l-2s„2)e2^^j. 


H H W+^ W-, 


2sn,^y 


H H Z^j^ Zy 


1 /II/ 

Or 2o 2y 


H W+^ Wp Z, 


1 ie^ fjixv 

2 c^y 


H W^ W-^ Z, 


1 ie^ 

2c^y 


W^+M Wp 


2 s„2y 


W+^ W-y Zp Zp 


2s„2y 


W+^ Wp Z, Zp 


2Cn,y 


W^ W-^ Z, Zp 


1 fii/ 

Of. y 


wt Wp z^ z. 


1 (l-2s^2)2g2 

2 Cw Sw 


Zf, Z,, Zi Zf 
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• integers: they consist of optional sign followed by one or more decimal digits, such as 

1 -1 123 -98765 
Integers can appear in Lagrangian terms, parameter definition and in other expressions. 

• Floating point numbers include optional sign, several decimal digits of mantissa with an 
embedded period (decimal point) with at least one digit before and after the period, and 
optional exponent. The exponent, if present, consist of letter E or e followed by an optional 
sign and one or more decimal digits. Valid examples of floating point numbers are 

1.0 -1.0 0.000511 5.11e-4 

Floating point numbers are used only as parameter values (coupling constants, particle 
masses etc). They can not be explicitly used in Lagrangian terms. 

• String constants may include arbitrary symbols. They are used as comments in parameter 
statements, full particle names in the declaration of a particle, etc. Examples from the 
previous section are 

electron 
photon 

If a string constant contains any character besides letters and digits or does not begin with 
a letter, it should be quoted. For example, the comments in QED and QCD input files 
(see previous section) contain blank spaces, so they are quoted: 

'elementary electric charge' 
'Strong coupling' 

2.2 Comments 

User can include comments into the LanHEP input file in two ways. First, symbol '%' denotes 
the comment till the end of current line. Second way allows one to comment any number of lines 
by putting a part of input file between ' /* ' (begin of comment) and ' */ ' (end of comment) 
symbols. 

2.3 Including files 

LanHEP allows the user to divide the input file into several files. To include the file file, the 
user should use the statement 

read Gle. 

The standard extension ' .mdl' of the file name may be omitted in this statement. 
Another way to include a file is provided by the use statement as 

use file. 

The use statement reads the file only once, next appearances of this statement with the same 
argument do nothing. This function prevents multiple reading of the same file. This form 
can be used mainly to include some standard modules, such as declaration of Standard Model 
particles to be used for writing some extensions of this model. 
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2.4 Conditional processing of the model 

Let us consider the LanHEP input files for the Standard Model with both the t'Hooft-Feynman 
and unitary gauges. It is clear that these input files differ by only a few lines — the declaration 
of gauge bosons and Higgs doublet. It is more convenient to have only one model definition 
file. In this case the conditional statements are necessary. LanHEP allows the user to define 
several keys, and use these keys to branch among several variants of the model. The keys have 
to be declared by the keys statement: 

keys namel^valuel, name2—value2, ... . 

Then one can use the conditional statements: 

do_if key==valuel. 

actionsl 
do_else_if key==value2. 

actions2 
do_else_if key==value3. 

do_else . 

default actions 
end_if . 

The statements do_else_if and do_else are not mandatory. The value of the key is a 
number or symbolic string. An example of using these statements in the Standard model may 
read: 

keys Gauge=unitary . 
do_if Gauge==Feynman . 

vector Z/Z: ('Z-boson' ,mass MZ = 91.187, width wZ = 2.502, gauge). 
do_else_if Gauge==unitary . 

vector Z/Z: ('Z-boson' ,mass MZ = 91.187, width wZ = 2.502). 
do_else . 

write ( 'Error : key Gauge must be either Feynman or unitary'), 
quit . 
end_if . 

Thus, to change the choice of gauge fixing in the generated Feynman rules it is enough to 
modify one word in the input file. There is another way, to set the key value from the command 
line at the launch of LanHEP: 

Ihep -key Gauge=Feynman filename. 

If the value of key is set from the command line at the program launch, the value for this 
key in the keys statement is ignored. 

There is also a short form of do_if statement: 

do_if (ifej==vaiue, (statementl), (statement2)). 

If the key is defined equal to value, then the statementl is processed. Otherwise the program 
executes statement2. The form with only one statementl is allowed. An example could be: 

do_if (electron_mass==on, (parameter Me=0 . 511e-3) , (let Me=0)). 
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3 Objects in the expressions for Lagrangian terms 

Each symbol which may appear in algebraic expressions (names of parameters, fields, etc) has 
a fixed order of indices and their types. If this object is used in any expression, one should 
write its indices in the same order as they were defined when the object was declared. 

Besides the types of indices corresponding to the color SU (3) group: color c3 (color triplet), 
color c3b (color antitriplet) and color c8 (color octet) described in the previous example, 
there are default types of indices for Lorentz group: vector, spinor and cspinor (antispinor). 
User can also declare new types of indices corresponding to the symmetries other than color 
SU (3) group. In this case any object (say, particle) may have indices related to this new group. 
This possibihty will be described in Section 9. 

If an index appears twice in some monomial of an expression, LanHEP assumes summation 
over this index. Types of such indices must allow the contraction, i.e. they should be one of 
the pairs: spinor and cspinor, two vector, color c3 and color c3b, two color c8. 

In general the following objects are available to appear in the expressions for a Lagrangian: 
integers and identifiers of parameters, particles, specials, let-substitutions and arrays. 

There are also predefined symbols i, denoting imaginary unit i {i^ — —1) and Sqrt2, which 
is a parameter with value \/2. 

3.1 Parameters 

Parameters are scalar objects {i.e. they have no indices). Parameters denote coupling con- 
stants, masses and widths of particles, etc. To introduce a new parameter one should use the 
parameter statement, which has the generic form 

parameter name=value : comment. 

• name is an identifier of newly created parameter. 

• value is an integer or fioating point number or an expression. One can use previously 
declared parameters and integers joined by standard arithmetical operators ' + ' , '-' , 
'*' , V , and ' ** ' (power) . 

• comment is an optional comment to clarify the meaning of parameter, it is used in Com- 
pHEP help windows. Comment has to be a string constant, so if it contains blank spaces 
or other special characters, it must be quoted (see Section 3). 

In the expression for parameters one can use functions sqrt , pow, sin, asin, cos, acos, 
tan, atan, atan2, f abs (same as in C programming language). One can add new functions 
through the use of the statement 

external_f unc (function, arity) . 

Arity is the number of arguments of the function. If this function is absent in the standard 
C hbrary, it should be written in C and added to the library file usrlib.a in the CompHEP 
working directory. 

3.2 Particles 

Particles are objects to denote physical particles. They may have indices. It is possible to use 
five statements to declare a new particle, at the same time statements define the corresponding 
Lorentz index(es): 
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scalar P/aP: {options) . 
spinor P/aP: (options) . 
vector P/aP: (.options) . 
spinor 3 P/aP: (options) . 
tensor P/aP: (options) . 

spinorS stands for spin 3/2, P and aP are identifiers of particle and antiparticle. In the case 
of truly neutral particles (when antiparticle is identical to the particle itself) one should use 
the form P/P with identical names for particle and antiparticle. 
It is possible to write only the particle name, e.g. 

scalar P : (options) . 

In this case the name of corresponding antiparticle is generated automatically. It satisfies the 
usual CompHEP convention, when the name of antiparticle differs from particle by altering the 
case of the first letter. So for electron name el automatically generated antiparticle name will 
be El. If the name contains symbol '+' it is replaced by '-' and vice versa. 

The option is comma-separated list of options for a declared particle, and it may include the 
following items: 

• the first element in this list must be the full name of the particle, (e.g. electron and 
photon in our example.) Full name is a string constant, so it should be quoted if it 
contains blanks, etc. 

• mass parani=value defines the mass of the particle. Here param is an identifier of a new 
parameter, which is used to denote the mass; value is its value, it has the same syntax as in 
the parameter statement, comment for this new parameter being generated automatically. 
If this option omitted, the mass is assumed to be zero. 

• width param^value declares the width of the particle. It has the same syntax as for mass 
option. 

• texname string declares the name to use in LaTeX output. 

• atexname string set LaTeX name for antiparticle, by default it is \bar{texname}. 

• itype is a type of index of some symmetry; one can use default index types for color SU (3) 
group (see QCD example in Section 2) . It is possible to use user-defined index types (see 
Section 9) and Lorentz group indices vector, spinor, cspinor. 

• left or right says that the massless spinor particle is an eigenstate of (1 — 75)/2 or 
(1 + 75)72 projectors, so this fermion is left-handed or right-handed. 

• gauge declares the vector particle as a gauge boson. This option generates the correspond- 
ing ghosts and goldstone bosons names for the named particle (see below) . 

When a particle name is used in any expression (in Lagrangian terms) , one should remember 
that the first index is either a vector or a spinor (if this particle is not a Lorentz scalar). Then 
the indices follow in the same order as index types in the options list. So, in the case of quark 
declaration (see the QCD example) the first index is spinor, and the second one is color triplet. 

There are several functions taking particle name as an argument which can be used in 
algebraic expressions. These functions are replaced with auxiliary particle names, which are 
generated automatically. 
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• Ghost field names in gauge tlieories are generated by tlie functions ghost (name) — > 
'name.c' and ccghost (name) — > 'name.C (see for instance Table 1). Here and below 
name is the name of the corresponding gauge boson. 

• Goldstone boson field name in the t'Hooft-Feynman gauge is generated by the function 
gsb (name) — > 'name.f. 

• The function ant i (name) generates antiparticle name for the particle name. 

• The name for a charge conjugate spinor particle iJj'^ = Cip'^ is generated by the function 
cc(name) 'name.c'. Charge conjugate fermion has the same indices types and or- 
dering, however index of spinor type is replaced by the index type cspinor and vice 
versa. 

• vev(expr) is used in the Lagrangian for vacuum expectation values. Function vev ensures 
that deriv*vev (expr) is zero. In other words, vev function forces LanHEP to treat expr 
as a scalar particle which will be replaced by expr in Feynman rules. 

3.3 Specials 

Besides parameters and particles other indexed such as 7-matrices, group structure constants, 
etc may appear in the Lagrangian terms. We refer such objects as specials. 
Predefined specials of the Lorentz group are: 

• delta(vector) can be used for g^^^ 

• gamma stands for the 7'^-matrices. It has three indices of spinor, cspinor and vector 
types. 

• gammaS denotes 75 matrix. It has two indices of spinor and cspinor types. 

• moment has one index of vector type. At the stage of Feynman rules generation this 
symbol is replaced by the particle moment. 

• deriv is replaced by —ip^, where p^ is the particle moment. It has one vector index. 
Note that deriv*A*B means {dA)B, not d{AB). 

• For the derivative of a product, derivp(ea;pr) can be used. derivp(A*B) is equivalent to 
deriv*A*B+A*deriv*B. 

• epsv is an antisymmetric tensor e'^^'^'^ with 4 tensor indices. 
Specials of the color SU{3) group are: 

• lambda denotes Gell-Mann A-matrices. It has three indices: color c3, color c3b and 
color c8. 

• f _SU3 is the SU{3) structure constant. It has three indices of color c8 type. 

• eps_c3 and eps_c3b are antisymmetric tensors. The first has 3 indices of the type color 
c3, and the second 3 color c3b indices. 

Note that for specials the order of indices types is fixed. 

Users can declare new specials with the help of a facility defined in Section 9 to introduce 
user-defined indices types. 
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3.4 Let-substitutions 

LanHEP allows the user to introduce new symbols and then substitute them in Lagrangian 
terms by some expressions. Substitution has the generic form 

let name=expr. 

where name is the identifier of newly defined object. The expression has the same structure as 
those in Lagrangian terms, however here expression may have free (non-contracted) indices. 
Typical example of using a substitution rule is a definition of the QED covariant derivative 

as 

let Deriv''mu=deriv~mu + i*ee*A"mu. 

corresponding to — 9^ + igeA^- 

More comphcated example is the declaration a'^^ = ^(7^7^ — 7*^7'*) /2 matrices: 

let sigma''a''b'"mu"nu = i*(gaiifflia'"a"c'"mu*gaiiima'"c'"b"nu 

- gaiifflia"a"c"nu*gainma''c"b''inu)/2. 

Note that the order of indices types of new symbol is fixed by the declaration. So, first two 
indices of sigma after this declaration are spinor and antispinor, third and fourth are vector 
indices. 

3.5 Arrays 

LanHEP allows to define components of indexed objects. In this case, contraction of indices 
will be performed as an explicit sum of products of the corresponding components. 
An object with explicit components has to be written as 

{exprl, expr2 exprN 

where expressions correspond to components. All indices of components (if present) have to 
be written for each component, and the index numbering components has to be written after 
closing curly bracket. Of course, all the components must have the same types of free (non- 
contracted) indices. 

Arrays are usually applied for the definition of multiplets and matrices corresponding to 
broken symmetries. 

Typical example of arrays usage is a declaration of electron-neutrino isospin doublet 11 (and 
antidoublet Ll) 

let ll~a~I = { nl"a, el~a}~I, Ll~a~I = { Nl"a, El~a}~I. 

Here we suppose that nl was declared as the spinor particle (neutrino), with the antiparticle 
name Nl. 

Note that the functions anti and cc can be applied also to the multiplets, so the construction 
let ll'^a^I = { nl''a, el"a}"I, Ll^a"^! = aiiti(ll) "a"^! . 
is possible. 

Matrices can be represented as arrays which have other arrays as components. However, it 
is more convenient to declare them with dummy indices, see Section 5.4 (the same is correct 
for multiplets also). 

It is possible also to use arrays directly in the Lagrangian, rather than only in the declaration 
of let-substitution. 

When LanHEP is launched, it has already declared some frequently used matrices. They 
are: 
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• deltaCAT) represents the unit matrix Sab of dimension N; 

• taul, tau2, tau3 are r- matrices ri, r2, T3: 



10 ' \ i or lo-i 



• tau is a vector r = (ti, T2, rs); 

• taup and taum are matrices — (r^ ± iT^)j\f2\ 

• taupm is a vector (r"*", r^, t~); 



eps is the antisymmetrical tensor e (s ^ = 1); 

Taul, Tau2, Tau3 are generators of SU{2) group adjoint representation (3-dimensional 
analog of r-matrices) T^,T'^,T^ with commutative relations [Tj,!}] = —ieijkTk'. 
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• Taup and Taum corresponds to = (T^ ± iT'^)/y/2; 

• Taupm is a vector f = T"). 

3.6 Two-component notation for fermions 

It is possible to write the Lagrangian using two-component notation for fermions. The con- 
nection between two-component and four-component notations is summarized by the following 
relations: 

'^-(i]- ^-(ir ^'-(if- 

If the user has declared a spinor particle p (with antiparticle P), the LanHEP notation for 
its components is: 

^ up(p) 

f] — > downCp) 

Tj — > up(cc(p)) or up(P) 

^ down(cc(p)) or down(P) 

To use the four-vector one should use the statement 

special sigma: (spinor2,spinor2, vector) . 

Note that the sigma object is not defined by default, thus the above statement is required. It is 
possible also to use another name instead of sigma (this object can be recognized by LanHEP 
by the types of its indices). 

LanHEP uses the following rules to convert the two-component fermions to four-component 
ones (we use Pr^l = (1 ± 75)/2): 
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V1V2 



^1-Pr^2 

riPRr2 

■ i'iPLr2 

■ ^tPR^2 



^ia^^2 = ^_17^^L^2 
fiicr^'r]2 = l/j^-f^PLl/^^ 



up(Pl)*up(p2) 
down (PI) *down (p2 ) 

up(pl)*up(p2) 
down(Pl)*down(P2) 

up(pl)*up(P2) 
down (pi ) *down (P2) 

up(Pl)*up(P2) 
down (pi) *down (p2 ) 

down ( P 1 ) * s i gma*up ( p2 ) 
down(pl)*sigma*up(P2) 



PI* (l-gammaS) /2*p2 

PI* (l+gammaS) /2*p2 

cc (pi) * (l-gammaS) /2*p2 

PI* (l+gammaS) /2*cc (P2) 

cc (pi) * ( l-gamma5) /2*cc (P2) 

cc (pi) * ( l+gammaS) /2*cc (P2) 

PI* (l-gammaS) /2*cc (P2) 

cc (pi) * (l+gammaS) /2*p2 

P 1 *gaiifflia* ( 1 -gammaS ) /2*p2 

cc (pl)*gamTna* ( l-gamma5) /2*cc (P2) 



4 Lagrangian expressions 

When all parameters and particles necessary for the introduction of a (physical) model are 
declared, one can enter Lagrangian terms with the help of the Iterm statement: 

Iterm expr. 

Elementary objects of expression are integers, identifiers of parameters, particles, specials, 
let-substitutions, and arrays. 

These elementary objects can be combined by the usual arithmetic operators as 

• exprl+expr2 (addition), 

• exprl-expr2 (subtraction), 

• exprl*expr2 (product), 

• exprl/expr2 (fraction; here expr2 must be a product of integers and parameters), 

• exprl**N {Nth power of exprl; iVmust be an integer). 

One can use brackets '(' and ')' to force the precedence of operators. Note, that indices can 
follow only elementary objects symbols, i.e. if Al and A2 were declared as two vector particles 
then valid expression for their sum is Al"mu+A2"mu, rather than (A1+A2) "mu. 



4.1 Where-substitutions 

More general form of expressions involves where-substitutions: 
expr where subst. 

In the simple form subst is naine=repl or several constructions of this type separated by 
comma In the form of such kind each instance of identifier name in expr is replaced by repl. 

Note that in contrast to let-substitutions, where-substitution does not create a new object. 
LanHEP simply replaces name by rcpl, and then processes the resulting expression. It means 
in particular that name can not have indices, although it can denote an object with indices: 

Iterm F**2 where F=deriv~mu*A~nu-deriv"nu*A~mu. 

is equivalent to 

Iterm (deriv"mu*A"nu-deriv"nu*A''mu) **2 . 



15 



The substitution rule introduced by the keyword where is active only within the current 
Iterm statement. 

More general form of where-substitution allows to use several name=repl substitution rules 

separated by semicolon In this case cxpr will be replaced by the sum of expressions; each 
term in this sum is produced by applying one of the substitution rules from semicolon-separated 
list to the expression expr. This form is useful for writing the Lagrangian where many particles 
have a similar interaction. 

For example, ifu,d,s,c,b,t are declared as quark names, their interaction with the gluon 
may read as 

Iterm gg*anti (psi) *gaiiima*lambda*G*psi where 

psi=u; psi=d; psi=s; psi=c; psi=b; psi=t. 

The equivalent form is 

Iterm gg*U*gamma*lambda*G*u + gg*D*gamma*lambda*G*d + 
gg*C*gamma*lainbda*G*c + gg*S*gamma*lambda*G*s + 
gg*T*gamma*lainbda*G*t + gg*B*gamma*lambda*G*b. 

Whcrc-substitution can also be used in let statement. In this case one should use brackets: 

let lsub= (expr yhere wsub=exprl). 

4.2 Adding Hermitian conjugate terms 

It is possible to generate hermitian conjugate terms automatically by putting the symbol 
AddHermConj to Iterm statement: 

Iterm expr + AddHermConj . 

Continuing the former example, one can write: 

Iterm a*H*H*h + b*H**3 + AddHermConj . 

Note, that the symbol AddHermConj adds the hermitian conjugate expressions to all terms 
in the Iterm statement. It means in particular that in the statement 

Iterm exprl + {expr2 + AddHermConj). 

the conjugate terms are added to both exprl and expr2. Thus, one should not place self- 
conjugate terms in the Iterm statement where AddHermConj present (or one should supply 
these terms with 1/2 factor). 

4.3 Using the superpotential formalism in the MSSM and its extensions 

In supersymmetric models (in particular, in the Minimal Supersymmetric Standard Model 
(MSSM) [5]) one makes use of the superpotential — a polynomial W depending on scalar fields 
Ai (superpotential also can be defined in terms of superfields, we do not consider this case). 
Then, there is the contribution to the Lagrangian: Yukawa terms in the form 



and F*Fi terms, where = dW/dA^ (for more details, see [5] and references therein). 

To use this formahsm in LanHEP, one should define first the multiplets of matter fields 
and then define the superpotential through the let-substitution statement. The example of the 
MSSM with a single generation may read: 
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keep_lets W. 

let W=eps*(mu*Hl*H2+inl*Hl*L*R+md*Hl*Q*D+mu*H2*Q*U) . 

Here symbols HI , H2, L, R, Q, U, D are defined somewhere else as doublets and singlets in 
terms of scalar particles. 

Note that before the definition of W this symbol should appear in the keep_lets statement. 
It is necessary to notify LanHEP that let-substitutions (multiplets) at the definition of W should 
not be expanded. Without this statement, the representation used by LanHEP for W will not 
contain symbols of multiplets but only the particles which were used at multiplets definition. 

Since W was declared in the keep_lets statement and contains the symbols of multiplets, 
one can evaluate the variational derivative of W by one or two multiplets, e.g. df (W,H1) or 
df (W,H1,L). Thus the Yukawa terms may be written: 

Iterm - df (W,Hl,H2)*fHl*fH2 - ... + AddHermConj . 

Here f HI , f H2 are fermionic partners of corresponding multiplets. 

To introduce the F*Fi terms one needs to declare the conjugate superpotential, e.g. Wc, and 
write: 

Iterm - df (W,Hl)*df (Wc.Hlc) - 

A better way is to use the function dfdfc(W,Hl) instead: 

Iterm - dfdfc(W,Hl) - 

The function evaluates the variational derivative, multiply it by the conjugate expression and 
returns the result. Moreover, it can introduce auxiliary fields to split vertices with 4 color 
particles (in the case of CompHEP output); see Section 10 for more details. 

4.4 Generation of counterterms for calculation at 1 loop 

Ultraviolet divergencies in renormalized quantum field theories are compensated by renormal- 
ization of wave functions (pi (f)i{l + SZi), masses Mj — >• Mi + 6 Mi, charges e.; — > Ci + 6ei, 
and may be some other Lagrangian paramters. Such transformation of the Lagrangian assumes 
changes in the Feynman rules. In particular vertices are changed due to the appearence of new 
terms — counterterms. It is possible to treat this transformation at 1-loop level by means of 
LanHEP package. Currently this feature is not supported by CompHEP, but it can be used for 
1-loop calculation by FeynArts [2] and FormCalc [3] (e.g. it is done in the SloopS [9] package). 

In 1-loop approximation renormalization parameters appear in expressions only at first 
power, i.e. SZiSZj = 0, where SZi means all renormalization parameters {5e, 5Mi, 5Zi, ...). 
This should be done with the help of the following statement 

infinitesimal d-Zl, d-Z2, d-Zn. 

It declares that a monomial should be omitted if it contains more than one of parameters d-Zl, 
d_Z2, d^Zn. 
The statement 

transform ohj -> expr. 

forces LanHEP to substitute symbol of a parameter or a particle, obj, by the expression expr 
when this object appears in the Lagrangian expressions of the Lagrangian. 

For example, if the electric charge e in QED is renormalized than the statement 

transform ee -> ee+dZe. 
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makes LanHEP substitute each occurence of ee symbol in further expressions by ee+dZe. 

The transform statement is very similar to the let statement and uses the same syntax 
for expression expr. It has however two advantages. First, one does not need to introduce new 

name for transformed parameter or field. Second, if one has the LanHEP file for some physical 
model, it is enough to add counterterms statements into the input file just after the statements 
declaring parameters and particles. 

An example for QED model with renormalization reads as: 

parameter ee = 0.3133: 'Electric charge', 
vector A/A: photon. 

spinor el/El : (electron, mass me=0. 000511). 
infinitesimal dZee, dZme, dZell, dZelr, dZA. 
transform ee -> ee*(l+dZee), me -> me+dZme, 
A -> A*(l+dZA), 

el -> el + (dZel/2*(l-gamma5)/2+dZer/2*(l+gamma5)/2)*el, 
El -> El + (dZel/2*(l+gamma5)/2+dZer/2*(l-gamma5)/2)*El. 
Iterm e*El* (i*gamma*deriv+gamma*A+me) *el . 

Sometimes it is necessary to block the shifts set by transform statement in some parts 
of the Lagrangian (depending on the specific renormalization scheme). This can be done by 
BlockTransf option: 

option BlockTransf =0 . 

statements... 

option BlockTransf— 1. 

Another option Inf iOrder can be used in the same way to set the maximal power of in- 
finitesimal parameters in the Lagrangian expressions. By default, it is set to 1. It is also 
possible to define high order renormalization parameters by applying **N in the declaration, 
like: 

infinitesimal d_Z12**2. 

If one use the two-component fermion notation, the shifts for the upper and the lower 
components of a 4-componcnt fcrmions must be defined. This is done automatically if the shift 
for 4-component fermion is defined as above for the electron, with left and right projectors. A 
special case concerns Majorana fermions, e.g. neutralino in MSSM. The shift should be defined 
here in a following way: 

spinor n/n: (neutralino, mass Mn) . 
infinitesimal dZn/dZnc. 

transform n -> n + (dZn/2*(l-gamma5)/2+dZnc/2*(l+gamma5)/2)*n. 

The second statement declares renormalization parameter dZn and its complex conjugate dZnc. 



4.5 Constructing the ghost Lagrangian 

The ghost Lagrangian can be constructed by means of the BRST formahsm. To use it, the user 
should first define the BRST transformations for fields (see Section 6.4). 
The gauge-fixing Lagrangian reads as: 

CGF^G-G+ + \\G'\' + \\a'\\ 
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where {i = ±, 7, Z) arc gauge-fixing functions. The ghost Lagrangian ensures the BRST 
invariance of the entire Lagrangian and can be written as ([8]) 

J^Gh — —^^BRStG^ + ^BRSTJ^Gh- 

where SbrstCgh is an overall function, which is BRST- invariant (at 1 loop we can set this to 
zero). 

So, for the photon and = {d ■ A), the LanHEP code for gauge-fixing and ghost terms 
reads as: 

let G_A = deriv*A. 

Iterm -G_A**2/2. 

Iterm - ' A . C ' *brst (G_A) . 

Here the brst function is used to get BRST-transformation of the specified expression. 

The inverse BRST transformation can also be used. One can declare the transformations 
for the fields by means of brsti_transf orm. The function brsti(expr) can be used in Iterm 
statements. 

5 Omitting indices 

Physicists usually do not write all possible indices in the Lagrangian terms. LanHEP also 
allows a user to omit indices. This feature can simplify the introduction of the expressions and 
makes them more readable. Compare two possible forms: 

Iterm El''a*gaifflna''a'"b''mu*A"mu*el''b. 

corresponding to geea{x)j^,,eb{x)Aij,{x), and 

Iterm El*gairana''mu*el*A''mu. 

corresponding to gee{x)'j^e{x)A^{x)). Furthermore, while physicists usually write vector indices 
explicitly in the formulas, in LanHEP vector indices also can be omitted: 

Iterm -i*ee*El*gamma*el*A. 

Generally speaking, when the user omits indices in the expressions, LanHEP faces two problems: 
which indices were omitted and how to contract indices. 

5.1 Restoring the omitted indices 

When the indexed object is declared the corresponding set of indices is assumed. Thus, if the 
quark q is declared as 

spinor q:('some quark', color c3) . 

its first index is spinor and the second one belongs to the color c3 type. If both indices are 
omitted in some expression, LanHEP generates them in the correspondence to order (spinor, 
color c3). However, if only one index is written, for example in the form q"a, LanHEP has 

to recognize whether the index a is of color c3 or of spinor types. 

To solve this problem LanHEP looks up the list of indices omitting order. By default this 
list is set to 

[spinor, color c3, color c8, vector] 
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The algorithm to restore omitted indices is the following. First, LanHEP assumes that 
the user has omitted indices which belong to the first type (and corresponding antitype) from 
this list. Continuing with our example with particle q one can see that since this particle is 
declared having one spinor index (the first type in the list) LanHEP checks whether the number 
of indices declared for this object without spinor index equals the number of indices written 
explicitly by user. In our example (when user has written q~a) this is true. In the following 
LanHEP concludes that the user omitted the spinor index and that explicitly written index is 
of color c3 type. 

In other cases, when the supposition fails if the user has omitted indices of the first type in 
the list of indices omitting order, LanHEP goes to the second step. It assumes that user has 
omitted indices of first two types from this fist. If this assumption also fails, LanHEP assumes 
that the user has omitted indices of first three types in the list and so on. At each step LanHEP 
subtracts the number of indices of these types assumed to be omitted from the full number of 
indices declared for the object, and checks whether this number of resting indices equals the 
number of explicitly written indices. If LanHEP fails when the list of indices omitting order is 
completed, error message is returned by the program. 

Note that if the user would like to omit indices of some type, he must omit all indices of 
this type (and antitype) as well as the indices of all types which precede in the list of indices 
omitting order. 

For example, if object Y is declared with one spinor, two vector and three color c8 indices, 
than 

• the form Y"a"b"c"d"e"f means that the user wrote all the indices explicitly; 

• the form Y"a"b"c"d"e means that the user omitted spinor index and wrote vector and 
color c8 ones; 

• the form Y~a~b means that the user omitted spinor and color c8 indices and wrote only 
two vector ones; 

• the form Y means that the user omitted all indices; 

• all other forms, involving different number of written indices, are incorrect. 

One could change the list of indices omitting order types by the statement SetDef Index. 
For example, for default setting it looks like 

SetDef Index (spinor , color c3, color c8, vector). 

Each argument in the list is a type of index. 

5.2 Contraction of restored indices 

A dummy index can be contracted only with some another dummy index. LanHEP expands 
the expression and restores indices in each monomial. LanHEP reads objects in the monomial 
from the left to the right and checks whether the restored indices are present. If such index 
appears LanHEP looks for the restored index of the appropriate type at the next objects. 
Note, that the program does not check whether the object with the first restored index has 
another restored index of the appropriate type. Thus, if F is declared as let-substitution for the 
strength tensor of electromagnetic field (with two vector indices) then expression F*F (as well 
as F**2) after processing dummy indices turns to the implied form F"mu"nu*F"mu~nu rather 
than F"mu''inu*F"nu''nu. 
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This algorithm makes the contractions to be sensitive on the order of objects in the monomial. 
Let us look again at the QED example. Expression El*gainma*A*el (as well as A*el*gainma*El) 
leads to correct result where the vector index of photon is contracted with the same index of 
7-matrix, spinor index of electron is contracted with antispinor index of 7-matrix and anti- 
spinor index of positron is contracted with spinor index of 7-matrix. However the expression 
El*el*A*gaiiima leads to the wrong form El"a*el~a*A~mu*gaiimia"b"c"mu, because the first 
antispinor index after the electron belongs to positron. Spinor indices of gamma stay free (non- 
contracted) since no more objects with appropriate indices (so, LanHEP will report an error 
since Lagrangian term is not a scalar). 

Note, that in the vertex with two 7-matrices the situation is more ambiguous. Let's look at 
the term corresponding to the electron anomalous magnetic moment e{x){'j^j^ — ■j^'y'^)e{x)Ffj^^. 
The correct LanHEP expression is 

el*(gamma"mu*gamma''nu - gamma''nu*gamma"mu)*El*F"mu"nu 

Here vector indices can't be omitted, since it lead to the contraction of vector indices of 7- 
matrices. One can see also that the form 

el*El*(gamma''mu*gamma"nu - gamma''nu*gamma~mu) *F"mu"nu 

will correspond to the expression e{x)e{x)Tr{'^^Y)PiJ.y Here LanHEP does see a scalar Lorentz- 
invariant expression in the Lagrangian term, so it has no reason to report an error. 

These examples mean that the user should clearly realize how the indices will be restored 
and contracted, or (s)he has to write all indices explicitly. 

5.3 Let-substitutions 

Another problem arises when the dummy indices stay free, this is the case for the let statement. 
LanHEP allows only two ways to avoid ambiguity in the order of indices types: either the 
user specifies all the indices at the name of new symbol and free indices in the corresponding 
expression, or he should omit all free indices. 

In the latter case the order of indices types is defined following the order of free dummy 
indices in the first monomial of the expression. For example if Al and A2 are vectors and cl 
and c2 are spinors, the statement 

let d=Al*cl+c2*A2. 

declares a new object d with two indices, the first is a vector index and the second is a spinor in- 
dex according to their order in the monomial Al*cl. Of course, each monomial in the expression 
must have the same type of free indices. 

5.4 Arrays 

The usage of arrays with dummy indices allows us to define matrices conveniently. For example, 
the declaration of r-matriccs (defined in Section 3.5) can be written as 

let taul = {{0, 1}, { 1, 0}}. 
let tau2 = {{0, i}, {-i, 0}}. 
let tau3 = {{1, 0}, { 0,-1}}. 

One can see that in such way of declaration a matrix is written " column by column" . 

The declaration of objects with three 'explicit' indices can be done using the objects already 
defined. For example, when r-matrices are defined as before, it is easy to define the vector 
T = (ti,T2,T3) as 
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let tau = {taul, tau2, tau3} . 

The object tau has three indices, first pair selects the element of the matrix, while the matrix 
itself is selected by the third index, i.e. tau~i~j~a corresponds to r,". 

On the other hand, the declaration of structure constants of a group is more complicated. 
Declaring such an object one should bear in mind that omitting indices implies that in a 

sequence of components the second index of an object is changed after the full cycle of the first 
index, the third index is changed after the full cycle of the second one, etc. For example, a 
declaration of the antisymmetrical tensor e"''^'^ can read as 

let eps = {{{0,0,0}, {0,0,-1}, {0,1,0}}, 
{{0,0,1}, {0,0,0}, {-1,0,0}}, 
{{0,-1,0}, {1,0,0}, {0,0,0}}}. 

One can easily see that the components are listed here in the following order: 

^111 ^211 ^311 ^121 ^221 ^321 ^131 ^233 ^333 

The declaration of more complex objects such as SU (3) structure constants can be made in the 
same way. 

6 Checking the correctness of the model 

6.1 Checking electric charge conservation 

LanHEP can check whether the introduced vertices satisfy electric charge conservation law. 
This option is available, if the user declares some parameter to denote elementary electric charge 
(say, ee in QED example), and then indicates, which particle is a photon by the statement 

SetEM (piioton, param). 

So, in example of Section 1 this statement could be 

SetEM(A,ee) . 

Electric charge of each particle is determined by analyzing its interaction with the photon. 
LanHEP checks whether the sum of electric charges of particles in each vertex equals zero. 
Command line option -v- charges allows to display the electric charge detected. 

6.2 Testing Hermitian conjugate terms 

LanHEP is able to check the correctness of the hermitian conjugate terms in the Lagrangian. 
To do this, user should use the statement: 

CheckHerm. 

For example, let us consider the following (physically meaningless) input file, where a charged 
scalar field is declared and cubic terms are introduced: 

scalar h/H. 
parameter a,b,c. 

Iterm a*(h*h*H+H*H*h)+b*h*h*H+c*H*H*h + h**3. 
CheckHerm. 

The output of LanHEP is the following (here 2 is the symmetry factor) : 
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CheckHerm: vertex (h, h, h) : conjugate (H, H, H) not found. 
CheckHerm: inconsistent conjugate vertices: 
(H, h, h) (H, H, h) 

2*a <-> 2*a 

2*b <-> (not found) 

(not found) <-> 2*c 

If the conjugate vertex is not found, a warning message is printed. If both vertices are present 
but they are inconsistent, more detailed output is provided. For each couple of the incorrect 
vertices LanHEP outputs 3 kinds of monomials: 1) those which are found in both vertices 
(these monomials are correctly conjugated), 2) the monomials found only in first vertex, and 
3) the monomials found only in the second vertex. 

6.3 Probing the kinetic and mass terms of the Lagrangian 

LanHEP allows user to examine whether the mass sector of the Lagrangian is consistent^. To 
do this, use the statement 

CheckMasses . 

This statement must be put after all Item statements of the input file. 

When the CheckMasses statement is used, LanHEP creates the file named masses . chk 
(in the directory where LanHEP was started). This file contains warning messages if some 
inconsistencies are found: 

• missing or incorrect kinetic term; 

• mass terms with mass different from the value specified at the particle declaration; 

• off-diagonal mass terms. 

Note that in more involved models hke the MSSM the masses could depend on other parame- 
ters and LanHEP is not able to prove that expressions in actual mass matrix and in parameters 

declared to be the masses of particles are identical. Moreover, it is often impossible to express 
the masses as formulas written in terms of independent parameters. For this reason LanHEP 
evaluates the expressions appearing in the mass sector numerically (based on the parameter 
values specified by the user) . 

It is typical for the MSSM that some fields are rotated by unitary matrices to diagonalize 
mass terms. In some cases, values of mixing matrices elements can not be expressed by formulas 
and need to be evaluated numerically. LanHEP can recognize the mixing matrix if the elements 
of this matrix were used in OrthMatrix statement. LanHEP restores and prints the mass 
original matrix before fields rotation by mixing matrix. 

The numerical check of correctness of the mass matrix diagonalization can be done if LanHEP 
has access to the routine which calculates the masses and the rotation matrix. To provide this 
access, one should should make the dynamic library, e.g. 

gcc -shared -o diag.so diag.c 
and then specify this library as the third argument of external_f unc statement: 

external_f unc (diag ,3,' . /diag . so ' ) . 
Note that './' is necessary unless you put your file into system shared library path. 

^ In CompHEP, the information about particles propagators is taken from the particle declaration, where particle mass and width 
(for Breit-Wigner's propagators) are provided. Thus, the user is not obliged to supply kinetic and mass terms in Iterm statements. 
Even if these terms are written, they do not affect the CompHEP output. 
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6.4 Checking BRST invariance 

LanHEP can check the BRST invariance (see the reviews in [6, 7]) of the Lagrangian. First, 
the user should declare the BRST transformations for the fields in the model by means of 
brst_transf orm statement: 

brst_transf orm field -> expression . 

For example, the BRST transformation for the photon in Standard Model SA^ — 8^0"^ + 
ie{W^c~ — W~c'^) can be prescribed by the statement: 

brst_transf orm A -> deriv*'A.c'+i*EE*('W+'*'W-.c'-'W-'*'W+.cO • 

The file 'sm_brst.mdl' in LanHEP distribution contains the code for gauge and Higgs fields 
transformation corresponding to the CompHEP implementation of the Standard model. 
Since the transformations are defined, the statement 

CheckBRST. 

enables the BRST transformation for Lagrangian terms (so it should be placed before the first 
Iterm statement). The output is CompHEP or LaTeX file. Certainly, if the Lagrangian is 
correct, the output files are empty. However some expressions identical to zero could be not 
simplified and remain in the output. 

6.5 Extracting vertices 

New statement allows to extract a class of vertices into separate file. This feature is useful in 
checking large models with thousands of vertices, such as the MSSM. The statement has the 
form 

SelectVertices(/iZe; vlist ) . 

Here file is a file name to output selected vertices, and vlist determines the class of vertices. 
The pattern vlist may have two possible formats. 
In the first format vlist is a list of particles: 

[PI, P2, P3, ... Pnl 

All vertices consisting only of the listed particles PI, P2,... are extracted. For example, the 
pattern 

[A, Z, 'W+', 'W-', el, El, nl, Nl] 

selects vertices of leptons of first generation interaction with gauge fields and the self-interaction 
of gauge fields (we used here particles notation of CompHEP). 
The second format reads 

\_P11/P12../Pln, P21/.../P2m, ... ] 

Here are several groups of particles, each group is joined by slash '/' symbol. Selected vertices 
have the number of legs equal to the number of groups in the list, each leg in the vertex 
corresponds to one group. For example, 

[el/El/nl/Nl, el/El/nl/Nl, A/Z/'W+'/'W-'] 
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selects vertices with two legs corresponding to leptons, third leg is gauge boson. Thus, this 
pattern extracts vertices of leptons interactions with gauge bosons. Note that the vertices of 
gauge bosons interaction are not selected in this case. If a group consists of one particle, it 
must be typed twice to enable second format (e.g. [Pl/Pl, P2/P2, P3/P3] for the vertex 
consisting of particles PI , P2, P3. 

It is possible to specify some options in SelectVertices statement. The WithAnti option 
adds antiparticle names for each particle in the list (in the first format) or in the each group (in 
the second format. The Delete option removes selected vertices from the LanHEP's internal 
vertices list, so these vertices are not written in output file IgrngA^.mdl (IgrngA'^.tex). This 
could be useful if SelectVertices statements are supposed to distribute all vertices which 
should be in the model, than some 'unexpected' vertices can be found in output IgrngA/^.mdl 
file. 

An example of lepton-gauge vertices selection could read as 
SelectVertices ( ' lept -gauge .mdl ' , [el/nl, el/nl, A/Z/'W+'] , WithAnti, Delete). 

7 Simplifying the expression for vertices 

7.1 Orthogonal matrices 

If some parameters appear to be the elements of the orthogonal matrix such as quark mixing 
Cabbibo-Kobayashi-Maskava matrix, one should declare them by the statement 

OrthMatrixC {{an, ai2, 013} , {021, 022, 023}, {031, 032, 033}} ). 

where aij denote the parameters. Such a declaration permits LanHEP to reduce expressions 
which contain these parameters by taking into account the properties of the orthogonal matrices. 
If the parameters are named as above, all , al2, ... a33 one can use the short form: 

OrthMatrix(a__ , 3) . 

Here two underline '_' characters stands for 2 digits, and the numer 3 sets the dimension of 
matrix. 

Note that this statement has no relation to the arrays; it just declares that these parameters 
a -ij satisfy the correspondent relations. If necessary, one can declare further a matrix with these 
parameters as components by means of the let statement. 

7.2 Reducing trigonometric expressions 

Some physical models, such as the MSSM and general Two Higgs Doublet Model, involve 
large expressions built up of trigonometric functions. In particlular, after the diagonalization 
of the Lagrangian in the Higgs sector, in the models mentioned above two angles, a and /3, 
are introduced, thus the Lagrangian may be written in LanHEP notation using the following 
definitions: 

parameter sa=0 . 5 : ' sinus alpha', 

ca=Sqrt (l-sa**2) :' cosine alpha', 
parameter sb=0.9: 'sinus beta', 

cb=Sqrt(l-sb**2) : 'cosine beta'. 

One can find in the output expressions like sa*cb+ca*sb = sin(Q; + /3) and much more com- 
plicated ones. To simplify the output and make it more readable the user can define new 
parameters, for example: 
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parameter sapb=sa*cb+ca*sb : 'sin(a+b) ' . 

In order to force LanHEP to substitute the expression sa*cb+ca*sb by the parameter sapb, 

one should use the statement 

SetAngle (sa*cb+ca*sb=sapb) . 

It is possible also to substitute an expressions by any polynomial involving paramters, for 
example, 

SetAngle ((sa*cb+ca*sb)**2+(sa*cb-ca*sb)**2=sapb**2+sainb**2) . 

where samb should be previously declared as a parameter. Note, that LanHEP expands both 
expressions before setting the substitution rule. The left-hand side expression can also contain 
symbols defined by let statement. In the Lagrangian expressions LanHEP expresses powers of 
cosines, ca**A'' and cb**A^ through sines using recursively the formula cos^ a = cos^^^(l — 
sin^ a) to combine all similar terms. The same procedure is performed for the left-hand side 
expression in SetAngle statement. 

If the substitution for one angle combination is defined, LanHEP offers to define the other 
ones for expressions consisting of parameters used in previous SetAngle statements. In our 
example, for all expressions consisting of parameters ca, cb, sa, sb, the following message 
is printed (e.g. for cos(a + f3)): 

Warning: undefined angle combination; use: 
SetAngle (ca*cb-sa*sb=aaOOO) . 

Here aaOOO is an automatically generated parameter. The message is printed only once for 
each expression. This feature is disabled by default; to enable producing these messages, one 
should issue the statement 

option Undef AngleComb=l . 



7.3 Heuristics for trigonometric expressions 

LanHEP can apply several heuristic algorithms to simplify these expressions. For each angle 
a, the user should declare parameters for sin a, cos a, sin2Q;, cos2q;, tana. Than one should 
use the angle statement: 

angle sin=pl, cos=p2, sin2=p3, cos2=p4, tan=p5, texname=i2a2i]c. 

Here pN — parameter identifiers, name — LaTeX name for angle, it is used to generate au- 
tomatically LaTeX names for trigonometric functions of this angle if these names are not set 
explicitly by SetTexName statement. This statement should immediately follow the declaration 
of the parameters for sin a and cos a. Only the sin and cos options are mandatory. Other 
parameters {i.e. sin2Q;, cos 2a, tana) should be declared if these parameters are defined before 
sin a and cos a. If the former parameters will be declared in terms of latter ones, they are 
recognized automatically and they need not appear in angle statements. 

For example, the declaration for trigonometric functions of P angle in MSSM may read: 

parameter tb=2 . 52 :' Tangent beta', 
parameter sb=tb/Sqrt (l+tb**2) : 'Sinus beta', 
parameter cb=Sqrt(l-sb**2) : 'Cosine beta'. 

angle sin=sb, cos=cb, tan=tb, texname=' \\beta' . 

parameter s2b=2*sb*cb: 'Sinus 2 beta', 
parameter c2b=cb**2-sb**2 : 'Cosine 2 beta'. 
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Here the parameters s2b and c2b are recognized automatically by LanHEP as sin 2a and cos 2a 
since they are declared in terms of sa, ca parameters. 

For a couple of angles ( a, /3 in this example) the user should declare the parameters for 
sin(Q; + (3), cos{a + (3), sm{a — (3), and cos(q; — (3) to allow using all implemented heuristics: 



These parameters arc recognized as trigonometric functions automatically, by analysis of the 
right-hand side expressions. 

It is possible to control the usage of heuristics by the statement: 

option SmartAngleComb=Ar. 

where AT is a number: 

heuristics are switched off; 

1 heuristics are switched on (by default) ; 

2 same as 1, prints the generated substitution rule if the simplified expressions consists of 
more than 3 monomials; 

3 same as 1, prints all generated substitution rule. 

4 same as 1, prints all generated substitution rule and some intermediate expressions (debug 



The substitution rules are printed as SetAngle statement and can be used for manual im- 
provement of expressions. 

7.4 Abbreviations 

The long expressions can be removed from the vertices and placed into the parameters table. 
This trick allows to optimize significantly the further work of matrix element computation both 
in memory consumption and computation speed. To generate these new auxiliary parameters 
(abbreviations) one should run LanHEP with the command line option -evl num. For the 
lengthly expression in a vertex (more then num monomials) an auxiliary variable is generated, 
and the expression is stored in the parameters table. nuni=2 is recommended. 

In case of FeynArts output, one should use -abbr command line option which produce 
optimal code for this program. 

8 Output 

Currently, LanHEP can generate Feynman rules in one of three formats: CompHEP, FeynArts, 



parameter sapb=sa*cb+ca*sb 

parameter samb=sa*cb-ca*sb 
parameter capb=ca*cb-sa*sb 
parameter camb=ca*cb+sa*sb 



'sin(a+b) ' . 
' sin(a-b) ' . 
' cos (a+b) ' . 
' cos (a-b) ' . 



mode) . 



LaTeX. 
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8.1 CompHEP 



By default, LanHEP produces 4 files with tables for independent parameters varsAT.mdl, de- 
rived ones funcAr.mdl, particles prtclsAT.mdl, and vertices IgrngAT.mdl, where AT is the num- 
ber specified in the model statement. To start working with CompHEP in the new model, one 
should copy these files into the models/ subdirectory of the CompHEP working directory. The 
number N should be next to the number of models already existing in this subdirectory. 

8.1.1 Modifying the particle table format 

New option allows to modify the format of the output particle table and to add new proprties 
(new columns in the table). The new format is set by means of prtcf ormat statement, for 
example, the default CompHEP table format can be set by: 

prtcf ormat fullname:' Full Name 
name : ' p ' , 

aname : ' ap ' , 

spin2, color, mass, width, aux, 
texname : ' latex P name ' , 
atexname : ' latex aP name ' . 

Each column of the table is described by the entry prop:title=value, where prop is the name 
of particle property, title is shown in the title line, and vaiue is default value used if the property 
of specific particle is not set. Important: the width of title fixes the width of the table column, 
so it should be wide enough to contain records for this property for any particle. If not specified, 
property name is used as title, and blank space for default value. 

There is a set of predefined properties (table columns) which exist in CompHEP by now 
plus the electric charge of a particle: 

• fullname 'Long' particle name. 

• name Particle name used in vertices table. 

• aname Antiparticle. 

• spin2 Twice spin, integer — 4. 

• mass Particle mass. 

• width Particle width. 

• echarge Particle electric charge, integer or ratio The value is generated automatically 
if the CheckEM(p/ioton, coupling) statement is used in the model file. 

• echargeS Three times electric charge, integer. 

• color Dimension of the color SU{3) group representation (one of 1,3,8). 

• aux Specify the particle as left or right fermion, boson in Feynman gauge (see CompHEP 
manual) . 

• texname LaTeX notation to be used for this particle. 

• atexname Same for antiparticle. If not set, barjtexname} is used. 

Besides these predefined properties, the user can introduce new ones. One can add, say, 
PDG particle number to the table: 
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prtcformat fullname:' Full Name 
name : ' p ' , 
aname : ' ap ' , 

spin2, color, mass, width, aux, 
pdg: 'PDG ID' , 

texname : ' latex P name ' , 
atexname : ' latex aP name ' . 

Then the pair prop value can be written in the particle declaration statement: 

scalar h:(higgs, mass Mh, pdg 123, width wh) . 

Another statement prtcprop can solve two problems. First, it can declare a property (or a 
few of them) which are not included in the particle table, but they still can be used in particle 
declarations for further use in the future. For example, one can declare pdg property: 

prtcprop pdg. 

and use it in the particle declarations even though CompHEP does not support this property 
now for future extensions. In this case, particle table format is not changed. Several properties 
can be listed, separated by comma. 

Another problem is to set a new property for all particles in the model. There is no need to 
modify all particles declarations in the model, one can use statement like this: 

prtcprop pdg:(h=123, g=124, f =125) . 

In this example h,g,f are some particles and the numbers are values of a property. Among the 
predefined properties, only the mass and width can be assigned (or changed) in this way. 

If some property is set to a number, a ratio, some particle or parameter name, then it can 
be used in the Lagrangian expressions as (prop particle), e.g. electric charge can be used to 
describe quark-photon interaction: 

Iterm (echarge q) *EE*Q*gamma*q*A . 

Of course one have to choose beween using the declared charge in expressions and automatic 
detection via SetEM statement. 

8.1.2 Lagrangian table setup 

The format and content of the Lagrangian (vertices) table can be tuned by several options set 
by option statement: 

option opname=value . 

and by some command line options. 

• chepCFWidth Width of Common Factor field, default is 50. 

• chepLPWidth Width of Lorentz Part field, default is 600. 

• RemDotWithFerm options (default value 1) tells LanHEP to replace gf^j, by {jixJu+JuJn)/'^ in 
the vertices with fermions (CompHEP requirement). If this replacement is not neccessary 
put this option to 0. 

• ReduceGammaS (default value 1) removes (1 ±75)72 operators in the vertices with lefthand 
or righthand fermions (e.g. neutrinos) replacing them by or 1. Put this option to to 
keep projector operators in such vertices. 
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• MultByl (default value for CompHEP output, 1 for LaTcX) In CompHEP tables, com- 
mon i in the vertices is not shown, so imaginary unit appears only in the vertices with 
pseudoscalars. Setting this option to 1 allows to restore it, making ComHEP vertices the 
same as in textbooks. 

• WriteAll (default 0) Setting this option to 1 makes LanHEP write into CompHEP table 

all vertices: 2-legs vertices, 1-legs ones (they can appear from incorrectly written Higgs 
potential) and also the vertices with more than 4 legs. The command line option -allvrt 

has the same effect. 

• MaxiLegs allows to limit the number of legs in the vertices produced when WriteAll is 
active. 

• WriteColors (default 0) Set this option to 1 to write color structure of vertices explicitely. 
Also the command line options -colors can be used. 

When the WriteColors option is set to 1, color matrices and dot products are written in 
the Lorentz Part, e.g. QCD plus quark-photon interactions produces the following vertices file: 



QCD 
Lagrangian 



PI 
G 



P2 

G 



IPS 

IG 



G.C 


iG.c 


IG 


Q 


Iq 


IG 


Q 


Iq 


lA 


G 


IG 


IG 



I 



P4 



> Factor <|> dLagrangian/ dA(pl) dA(p2) dA(p3) 

gg |m2.p3*ml.m3*F(cl jC2,c3) 

-ml . p3*m2 . m3*F (cl , c2 , c3) 
+m3.pl*ml.m2*F(cl,c2,c3) 
-m2 . pl*inl . m3*F (cl , c2 , c3) 
-m3.p2*inl.m2*F(cl,c2,c3) 
+ml . p2*m2 . in3*F (c 1 , c2 , c3) 

-gg |m3.p2*F(cl,c2,c3) 

gg |L(cl,c2,c3)*G(m3) 

ee/3 |cl.c2*G(m3) 

gg'^2 I ml . m3*m2 . m4*F ( c 1 , c2 , cO) *F ( c3 , c4 , cO) 

-ml . m4*m2 . m3*F (cl , c2 , cO) *F (c3 , c4 , cO) 
+ml . m2*m3 . m4*F (c 1 , c3 , cO) *F (c2 , c4 , cO) 
-ml .m4*m2 .m3*F (cl , c3 , cO) *F (c2 , c4 , cO) 
+ml.m2*m3.m4*F(cl,c4,c0)*F(c2,c3,c0) 
-ml . m3*m2 . m4*F (cl , c4 , cO) *F (c2 , c3 , cO) 



Here q/Q is a quark, G - gluon , G.c - gloun ghost field, and A - photon. F stand 

for antisymmetric structure constants (D for symmetric ones), L - Gell-Mann matrices. The 
Lorentz part is shown one monomial per line just to make the expression fit into the page, 
usually there is one line per vertex (option chepBreakLines set to 1 to made this effect). 



8.2 FeynArts format 

Feynman rules in FeynArts format are generated if the command line option -f eynarts is used 
at LanHEP start. In this case 4 files are generated. First, modeliV.gen contains 'generic 
model', in particular the Lorentz structures which appears in the model. N is the num- 
ber defined in the model statement in the LanHEP input file. Note that default generic 
model FeynArts/Models/Lorentz . gen can not be used with LanHEP output. Second file, 
modeliV.mod 'classes model' defines interaction vertices. Two other files are used to compile 
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FORTRAN program from FormCalc output. modelAT.h contains declarations of parameters 
and common blocks, mdl_iniiV.F serve to initialize parameters. 

To use the model generated by LanHEP one should have in Mathematica input file the 
statement like: 

SetOptions [InsertFields,Model->modeliV, GenericModel->modeliV] 

The process can be set, for example for electron-positron annihilation into two photons : 

process = {prt ["el"] ,prt ["El"] } -> {prt ["A"] ,prt ["A"] } 

Here particle names defined in the LanHEP model can be used. 

When FORTRAN code is generated by FormCalc, in the output replace the file model. h 
by modeliV.h, and in the file run.F find #include "sm_ini.F" and replace sin_ini.F by 
mdl_iniiV.F 

Note that FormCalc has some model specific definitions, which can be seen in the end of 
FormCalc. m file. They can change some parameters of the model generated by LanHEP. The 
safest way is to create model satisfying these definitions. In particular, for each particle mass 
(say MH) one should declare (in the LanHEP input file) its square like 

parameter MH2=MH**2. 

The command line option -abbr allows to create abbreviations for the interaction vertices, 
placing the expressions from vertices into parameters definition and making the model more 
compact and increasing the speed of calculation of the matrix element. 

If one has in mind 1-loop computations, rcnormalization parameters should be defined. One 
can put the expression to be used by FormCalc in the infinitesimal statement, e.g. for 
electric charge, weak angle and tadpole counterterm one can write: 

infinitesimal dEE= -(dZAA - SW/CW*dZZA)/2, 

dSW= CW"2/SW/2*(dMZsq/MZ"2 - dMWsq/MW~2) , 
dTad='-ReTilde [Self Energy [prt ["H"] ->{}, -1]]'. 

The expression should be in quotes " if it does not satisfy LanHEP syntax ('some [thing]' and 
'{}' do not). If the expression is not provided for a parameter, some default value is used as in 
the standard FcynArts SM.mod file, thus for the Standard Model only these 3 parameters have 
to be set cxplicitcly. Note that the old sign definition (as in FormCalc 4.1) is used for vector 
self-energies in the default scheme. 

This output mode was tested with FeynArts 3.2 and FormCalc 4.1. 

8.3 LaTeX output 

LanHEP generates LaTeX output instead of CompHEP model files if the user set -tex in the 
command line to start LanHEP. Three files are produced: 'varsJV.tex', 'prtclsiV.tex' and 
'lgrngiV.tex\ The first file contains names of parameter used in physical model and their 
values. The second file describes the particles, together with propagators derived from the 
vertices. 

The last file lists introduced vertices. LanHEP uses Greek letters u, p... for vector indices, 
letters a, b, c... for spinor ones and p, q, r... for color indices (and for indices of other groups, if 

they were defined). 

It is possible to inscribe names for particles and parameters to use them in LaTeX output. 
It can be done by the statement 
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SetTexNcone ( [ idcnt=tcxnaine, ... ]). 

Here ident is an identifier of particle or parameter, and texname is string constant containing 
LaTeX command. Note, that for introducing backslash '\' in quoted string constant one should 
type it twice: '\\'. 

For example, if one has declared neutrino with name nl (and name for antineutrino Nl) than 
the statement 

SetTexName([nl='\\nu~e' , Nl=' \\bar{\\nu}"e'] ) . 

makes LanHEP to use symbols and for neutrino and antineutrino in LaTeX tables. 
The vertices table can be tuned by the following command line options: 

-texLines num Set number of lines in LaTeX tables to num. After the specified number 
of lines, LanHEP continues writing current table on the next page of LaTeX the output. 
Default value is 40. 

-texLineLength num Controls width of the Lagrangian table. Default value is 35, user 
can vary table width by changing this parameter. 

-nocdot removes \cdot commands between the parameters symbols in the LaTeX output. 
This option is useful when all parameters have prescribed LaTeX names. 

-f rc If -tex option is set, forces LanHEP to split 4-fermion and 4-color vertices just as it 
is made for CompHEP files. 

9 Declaration of new index types and indexed objects 
9.1 Decleiring new groups 

Index type is defined by two keywords^: group name and representation name. Thus, color 
triplet index type color c3 has group name color and representation name c3. 
LanHEP allows user to introduce new group names by the group statement: 

group gname. 

Here gname is a string constant, which becomes the name of newly declared group. 
Representation names for each group name must be declared by the statement 

repres gname: (.rlist) 

where rlist is a comma-separated list of representation names declaration for the already de- 
clared group name gname. Each such declaration has the form either rnamc or rname/ crname. 
In the first case index which belongs to the gname rname type can be contracted with another 
index of the same type; in the second case index of gname rname type can be contracted only 
with an index of gname crname type. 

For example, definition for color SU{3) group with fundamental, conjugate fundamental and 
adjoint representations looks like: 

group color:SU(3). 
repres color: (c3/c3b,c8) . 

So, three index types can be used: color c3, color c3b, color c8. The contraction of 
these indices is allowed by pairs (color c3, color c3b) and (color c8, color c8) indices. 

^The exception is Lorentz group, corresponding indices types are defined by a single keyword. 
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9.2 Declaring new specials 

Specials with indices of user-defined types can be declared by means of special statement: 
special name: Qlist) . 

Here name is the name of new symbol, and Hist is a comma-separated list of indices types. For 
example, Gell-Mann matrices can be defined as (although color group and its indices types are 
already defined): 

special lambda: (color c3, color c3b, color c8) . 

To define Dirac's 7-matrices one can use the command 

special gamma: (spinor, cspinor, vector). 

9.3 Arrays 

Array, i.e. the object with explicit components, can also have the user-defined type of index. 
In this case generic form of such object is 

{ exprl, expr2 ... ,exprN ; itype } 

where N expressions exprl ... exprN of N components are separated by comma, and itype is 
an optional index type. If itype is omitted LanHEP uses default group name wild and index 
type wild N, where AT is a number of components in the array. 

10 Splitting the vertices with 4 colored particles 

The CompHEP Lagrangian tables do not describe explicitly the color structure of a vertex. 
If color particles are present in the vertex, the following implicit contractions are assumed 
(supposing p, q, r are color indices of particles in the vertex): 

• Spq for two color particles A^, A'^; 

• Xpg for three particles, which are color triplet, antitriplet and octet; 

• fP'i'^ for three color octets. 

Other color structures are forbidden in CompHEP. 

So, to introduce the 4-gluon vertex /^^''G'^G'^/^**G*G^ one should split this 4-legs vertex 
into 3-legs vertices JP'^^'GIGIXp^: 

G-. .-G G :----G G :.. ..-G 

'\-^-\ + x\ + x\ ">:;■ 

G-" -G G ■■■■■G G "■■■G 



Here the field X^^^^ is a Lorenz tensor and color octet, and this field has constant propagator. 
If gluon name in CompHEP is ' G ' , the name ' G . t ' is used for this tensor particle; its indices 
denoted as 'm_' and 'M_' [' is the number of the particle in table item). 

The described transformation is performed by LanHEP automatically and transparently for 
the user. Each vertex containing 4 color particles is split to 2 vertices which are joined by an 
automatically generated auxiliary field. 
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10.1 Optimization in case of supersymmetry 

The same technique is apphed in the MSSM where more vertices with 4 color particles appear: 
vertices with 2 gluons and 2 squarks and vertices with 4 squarks. However, the large amount 
of vertices with 4 squarks requires many auxiliary fields, which can easily break CompHEP 
limitations on the particles number. It is possible however to reduce significantly the number 
of vertices and auxiliary fields if one introduce auxiliary fields at the level of multiplets. 

The vertices with 4 squarks come from DD and F*F terms. For example, there is the term 
^DqDq in the Lagrangian, 

where Q, D, U are squarks multiplets, and A is Gcll-Mann matrix. Instead of evaluating this 
expression and that splitting all vertices independently, one can introduce one color octet 
auxiliary field and write this Lagrangian term as Dq^'^. 

Other DD terms contain both color and colorless particles. Thus, the term D\D\ with 

D\ = gi{Q*TQ + L*TL + H*T'Hi + H;TH2), 

can be represented as 

gi{Q*TQ)C + gf{Q*T'Q){L*TL + H*THi + H*T'H2) + gf{L*TL + H*T'Hi + H*TH2f 

where ^* is the triplet of auxiliary fields. This terms can also be written in the another form: 

gi{Q*T'Q + L*T'L)C + gl{Q*T'Q + L*TL){H*THi + H;T'H2) + gl{H*THi + H;TH2f, 

where all vertices with 4 scalars (except vertices with Higgs particles) are splitted. Although 
the latter splitting is not mandatory, it can reduce significantly the amount of vertices. 

The similar technique is applicable to the F*F terms, with the transformation F*Fi — > 
F*^ + F,ii. 

Thus, we distinguish two types of vertices splitting: splitting at multiplet level and splitting 
at vertices level. Note that splitting the vertices with two gluons and two squarks must be done 
at vertices level after combining the similar terms, otherwise they would contain the elements 
of squark mixing matrices. 

The vertices splitting at multiplet level is implemented in LanHEP mainly for MSSM needs. 
The first case refers to DD terms. The user should declare several let-substitutions and then 
put in Iterm statement the squared sum: 

let al=g*Q*tau*q/2, 

a2=g*L*tau*l/2, 

a3=g*Hl*tau*hl/2, 

a4=g*H2*tau*h2/2 . 
Iterm - ( al + a2 + a3 + a4 ) ** 2 / 2. 

In this case LanHEP looks for the square of the sum of several let-substitution symbols, each 
containing two color or merely scalar particles. If such an expression is found, it is replaced as 
in the previous formulas. 

The vertices splitting in F*F terms is performed by dfdf c function (see previous section). 
After taking the variational derivative the monomials with two color or scalar particles (except 
Higgs ones) are multiplied by auxiliary fields, thus mediating the vertices with 4 color (scalar) 
particles. 

The multiplet level vertices splitting is controlled by the statement 
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option SplitColl=iV. 
where iV is a number: 
-1 remove all vertices with 4 color particles from Lagrangian; 

turn off multiplet level vertices splitting; 

1 allows vertices splitting with 4 color multiplets; 

2 allows vertices splitting with any 4 scalar multiplets except Higgs ones (more generally, 

any multiplets containing vev's). 

The value of this option can be set to different values before executing different Iterm state- 
ments. 

The vertices level splitting is performed after combining similar terms of the Lagrangian. 
This sphtting can be controlled by the statement 

option SplitCol2=iV. 

where iV is a number: 

disable vertex level splitting; 

1 enable vertex level splitting (only for vertices with 4 color particles). 

For CompHEP output, the default value is 2 for SplitColl and 1 for SplitCol2. For 
FeynArts and LaTeX output, default value is for both options. 

10.2 Splitting vertices manually 

Sometimes it is necessary to split the vertex by hands. Usually this happens if the automatically 
generated vertices appears to be non-hermitian. Then one should write the Lagrangian using 
this auxiliary fields explicitely. An auxiliary particle can be generated by means of AuxPrt (.spin, 
color, charged) function. Here spin is one of 0, 1, or 2, color is one of 1, 3, 8, and charged is 
for truly neutral particle, and 1 otherwise. For charged particle one should write like 

let a=AuxPrt(0,l,l) , A=anti(a). 

Note that ant i (AuxPrt ( . . . )) will generate another auxiliary particle and thus is wrong. 

11 Installation and running LanHEP 

To install LanHEP on your computer, you should get the archive file from the WWW-page 
http://theory.sinp.msu.ru/~semenov/lanhep.html. Unpack the archive: 

gunzip lhep300.tar.gz 
tar xf IhepSOO.tar 

The archive contains the directory lanhep with C source files. If you do not have gcc compiler, 
you should tune the makefile, change 'gcc' to the compiler which you have. To create the 
executable file (called Ihep) , go to lanhep directory and type 

make 

When LanHEP is complied, remove the source files by typing 
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make clean 



The archive also contains the subdirectory mdl with startup file and examples for several phys- 
ical models. Add the directory containing LanHEP to your PATH environment variable. Then 
LanHEP can be started from any other directory, it can find automatically files in the mdl 
directory. 

As mentioned above, LanHEP can read the model description from the input file prepared 
by the user. To start LanHEP write the command 

Ihep filenaine options 

where the possible options are described in the next section. If the fiiename is omitted, LanHEP 
prints it's prompt and waits for the keyboard input. In the last case, user's input is copied into 
the file Ihep . log and can be inspected in the following. To finish the work with LanHEP, 
type 'quit. ' or simply press ~D (or "Z at MS DOS computers). 

Possible options, which can be used in the command line to start LanHEP are: 

-V Verbose output: to provide information about processing the model, -vv and -vvv 
provide more information. 

-OutDir directory Set the directory where output files will be placed. 

-InDir directory Set the default directory where to search files, which included by read 
and use statements. 

-c3 produces the output files for CompHEP version 3x.xx, which has slightly different 
format than 4x.xx. 
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